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I.  COMPUTATIONAL  STUDIES  OF  LONG-RANGE 
HIGH-FREQUENCY  IONOSPHERIC  DUCTING 

1. 1 Introduction 

Ray-tracing  calculations  for  HF  propagation  studies  using  various  three 
dimensional  ionospheric  models  for  rays  originating  from  the  earth  has 
shown  the  existence  of  a small  number  of  long-range  ducted  modes  which 
are  a strong  function  of  the  time  of  day  and  season.  A ducted  mode  is  one 
in  which  HF  electromagnetic  waves  traverse  long  distances  by  successive 
reflections  from  various  height  regions  of  the  ionosphere  without  reaching 
the  ground.  For  some  conditions  they  do  not  exist  at  all  while  for  others 
they  appear  with  a small  efficiency  relative  to  all  other  modes.  The  search 
for  these  ducted  rays  in  a given  ionosphere  by  searching  through  many  values 
of  the  various  circuit  parameters  does  not  yield  the  kind  of  results  which 
would  be  useful  in  predicting  the  general  conditions  necessary  to  support 
these  modes. 

In  order  to  have  greater  flexibility  in  exploring  ionospheric  conditions 
we  have  chosen  to  work  with  a two  dimensional  ionospheric  model  where  the 
individual  ionospheric  profiles  could  be  constructed  to  reflect  more  arbi- 
trarily specified  horizontal  and  vertical  electron  density  gradients.  For  the 
assessment  of  the  effects  of  ionization  gradients  on  ray  injection  and  ducting 
at  HF  frequencies  we  assume  that  both  the  effects  of  the  earths  magnetic 
field  and  the  interaction  between  the  electrons  and  neutral  atoms  can  be 
neglected.  We  also  assume  that  lateral  deviations  do  to  transverse  electron 
density  gradients  are  small. 

Some  results  on  the  dependence  of  the  ducting  efficiency  upon  various 
gradients  have  been  presented  in  an  earlier  paper.  ^ In  this  report  we 
expand  on  these  results  and  extend  the  discussion  to  include  some  additional 
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work  concerning  the  stability  of  these  modes  to  both  random  and  determi- 
nistic fluctuations  of  the  ionosphere.  We  also  include  a discussion  of  these 
long-range  modes  for  an  elevated  source. 


I.  2 Ray  Equations 

The  six  canonical  equations  in  spherical  coordinates  for  a ray  propagating 
in  an  ionosphere  without  magnetic  field  can  be  reduced  to  two  equations  if  the 
medium  is  confined  to  vary  only  in  two  dimensions.  The  propagation  geometry 
corresponds  then  to  a spherical  system  where  the  ray  is  confined  to  a plane. 
These  two  differential  equations  have  the  form 


dx  a 

do  _ I I tan_q  Sy.  (I_2) 

dx  ~ vi  5h  a H 3x 

where  a is  the  radius  of  the  earth,  h is  the  height,  r is  the  distance  of  a point 
on  the  ray  from  the  center  of  the  earth,  x is  the  angular  distance  measured 
from  the  starting  point  of  the  ray  at  the  earth's  surface,  o is  the  angle  between 
the  ray  direction  and  the  tangent  to  the  earth.  For  the  ionosphere  without 
magnetic  field  the  refractive  index  \i  is  related  to  the  plasma  frequency  fN  and 
the  propagation  frequency  f by 


2 *N  (x'h) 

H (x.h)  = 1 - ^ 


(1-3) 


Equations,  resulting  from  Eqs.  (1)  and  (2)  after  a change  of  variables, 
are  integrated  analytically  for  \i  = 1.  For  |J<1  they  are  integrated  numerically 
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using  a standard  fourth- order  Runge-Kutta  technique  for  two  simultaneous 

2 

differential  equations.  For  ray  computations,  a 10-km  stepping  interval 
in  ground  range  is  used. 


I*  3 Two-Dimensional  Ionospheric  Model 

For  mathematical  and  numerical  convenience  and  without  implying  any 
unusually  superior  relationship  to  physical  reality,  electron-density  height 
profiles  are  formed  using  three  segments  of  sine- square  functions  of  the  type 

N . N,  + (N2-Nl)Sin2  (1-4, 

specifying  the  relationship  between  electron  density  (N)  and  height  (h).  These 
profiles  are  changeable  through  the  choice  of  critical- frequency  parameters 
without  causing  discontinuities  in  vertical  ionization  gradients.  Ionospheric 
profiles  are  parametrically  specified  in  selected  increments  of  distance  along 
the  surface  of  a curved  earth  while  intermediate  values  of  electron  density 
are  obtained  by  interpolation  using  second-order  polynomials.  With  the  choice 
of  these  parameters,  vertical  and  horizontal  ionization  gradients  may  be  in- 
troduced into  an  ionospheric  cross-section  to  which  radio  rays,  specified  by 
'operating'  frequency  and  initial  elevation  angle,  are  confined. 

A profile  at  any  given  point  along  the  great  circle  path  is  entirely  described 
by  specifying  the  following  parameters: 

hQ  t*le  bottom  of  the  ionospheric  E layer, 

and  height  and  corresponding  plasma  frequency  of  the  E layer 

maximum, 

hFl  "d  £F1  the  height  and  corresponding  plasma  frequency  of  the  FI  layer 
maximum, 

hF2  and  fF2  4116  height  and  corresponding  plasma  frequency  of  the  F2  layer 
maximum. 
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Any  of  the  individual  layers  can  be  omitted  from  the  profile  by  setting 

the  parameters  of  the  layer  equal  to  the  ones  of  the  layer  below  it.  For 

example,  if  h = h and  f = f„  there  will  be  no  FI  layer.  One  feature 

of  the  model  is  that  the  vertical  gradient  of  the  electron  density  is  zero  at 

each  of  the  layer  heights?  The  interpolation  between  any  of  the  given  profiles 

is  such  that  N(x,h),  and  are  continuous  throughout  the  entire  region. 

Oil  ox 

For  most  of  the  ionospheres  used  in  this  study  the  profiles  were  specified 

at  2000  km  intervals  in  ground  range.  This  spacing  was  adequate  to  describe 

a variety  of  horizontal  and  vertical  electron  density  gradients.  In  order  to 

determine  the  capability  of  an  individual  ionosphere  to  support  long-range 

modes  we  extend  the  electron  density  model  over  a 10,000  km  distance  along 

the  surface  of  the  earth.  For  approximately  the  first  5000  km  the  horizontal 

3N 

ionozation  gradients  were  made  to  be  negative  and  thereafter  positive. 

This  allowed  for  a comparison  of  the  ground  range  for  various  ray  parameters 
since  those  rays  which  were  first  injected  into  an  ionospheric  duct  were 
eventually  reflected  to  the  ground  by  the  symmetric  gradients. 

Some  features  of  the  profiles  are  common  to  all  the  models.  The  iono- 
spheric E layer  starts  at  hQ  = 90  km  and  has  a maximum  ionization  at 

h =110  km.  While  the  value  of  this  maximum  varies  from  model  to  model 
E 

it  is  always  taken  to  be  constant  over  the  ray  path  so  that  there  is  no 

horizontal  electron  density  gradient  below  110  km.  When  an  FI  layer  is  used 

the  plasma  frequency  is  specified  at  h^  = 160  km  and  the  F2-layer  maximum 

is  located  at  h = 300  km.  The  choice  of  constant  heights  for  h and  h 

F2  FI  F2 

in  no  way  restricts  the  kinds  of  ionization  gradients  that  we  can  study  since 

dN 

the  net  — — at  any  point  in  space  can  still  be  controlled  by  the  choice  of  values 

OX 

of  the  plasma  frequencies  along  the  path. 

Parameters  for  many  of  the  ionospheres  that  were  studied  are  shown  in 
TABLE  I.  1.  The  labels  indicate  some  of  the  features  of  the  ionosphere.  The 
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Parameters  for  Two  Dimensional  Ionospheric  Models 
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number  refers  to  the  horizontal  gradient  due  to  the  change  in  foF2  as  a 

function  of  ground  range  with  1 being  the  smallest  and  3 the  largest  horizontal 

gradient.  The  letters  specify  one  of  four  possible  values  of  f^,  and  f^,  the 

first  refers  to  f and  the  second  to  f . Some  typical  profiles  are  shown  in 
1 E 

Figure  1-1  where  foF2  = 8.  0 MHz  in  all  cases.  The  electron  density  in  three 
of  the  profiles.  A'  , D1  , D are  monotonic  functions  of  the  height  and  A,  B,  and 
C have  ionization  valleys  at  160  km.  The  values  the  horizontal  ionization 
gradients  at  several  heights  for  the  three  models  is  shown  in  TABLE  1-2.  A 
large  variety  of  ionospheric  conditions  can  be  studied  by  combining  the 
separate  models  for  the  horizontal  and  vertical  electron  density  gradients. 


I.  4 Ducting  Efficiency  for  Various  Ionospheric  Models 

In  order  to  determine  the  relative  efficiency  for  ducting  for  each  of  the 
ionospheres  we  trace  an  array  of  rays  in  0.  5 MHz  steps  in  the  frequency 
domain  between  6 and  25  MHz  and  in  0.  5°  steps  in  ground  elevation  between 
0°  and  11°.  In  soma  cases  we  reduced  the  size  of  the  frequency  interval  in 
an  attempt  to  determine  the  sensitivity  of  finding  ducted  rays  as  a function 
of  the  sampling  size.  For  all  of  the  profiles  shown  in  Figure  1-1  no  long 
range  ducted  rays  were  found  when  horizontal  gradients  were  omitted. 

This  was  the  case  even  when  the  frequency  sampling  interval  was  reduced 
to  0.  1 MHz.  Each  ray  is  traced  until  one  of  the  following  conditions  is 
fulfilled:  the  ray  touches  the  ground,  the  local  elevation  angle  is  0°  in  a 
non-ionized  portion  of  the  path,  the  ray  penetrates  at  h = 400  km.  , or  when 
a ducted  ray  reaches  the  limit  in  range  over  which  the  ionosphere  is 
specified. 

For  model  3d'  22  rays  were  found  to  be  trapped,  reaching  ground 
ranges  in  excess  of  8000  km  before  touching  the  earth.  The  behavior  with 
distance  and  height  of  some  of  these  rays  is  illustrated  in  Figure  1-2, 
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TABLE  1-2 


Horizontal  Electron  Density  Gradients 


Ionospheric 

Models 

dN 

— at  x = 
dx 

2000  km.  (el.  /cc.  / m.  ) 

Height 

300  km 

200  km 

185  km 

1 (A,  B,  C,  and  D) 

-0.  05 

-0.0093 

-0. 0038 

2 (A,  B,  C,  and  D) 

-0.099 

-0.019 

-0. 0076 

3 (A,  B,  C,  and  D) 

-0. 140 

-0. 026 

-0.011 

together  with  the  iso-ionic  contours  of  the  modeled  ionosphere.  The  contours 

_5 

are  identified  by  number  of  el/cc  x 10  . Four  of  the  22  ducted  rays  are 

shown  for  elevation  angles  of  0°  and  10°.  For  an  initial  elevation  angle  of  10°, 
for  example,  rays  computed  for  13.5  and  14.5  MHz  return  to  the  ground  at 
relatively  short  distances,  while  a ray  for  14.0  MHz  appears  trapped  in  an 
ionospheric  channel.  For  zero- elevation  angle,  two  rays  at  19-5  and  20.0 
MHz  reach  large  distances,  while  rays  of  19  and  20.5  MHz  reach  a short  and 
moderate  distance  respectively.  The  nature  of  the  paths  for  the  two  elevation 
angles  is  quite  different.  The  10°  ray  has  five  oscillations  and  is  constrained 
to  a 40  km  height  interval  while  the  0°  ray  has  three  oscillations  over  a 90  km 
height  interval.  The  ducted  rays  occur  for  frequencies  just  above  those  for 
penetration  of  the  E layer  maximum.  A detailed  plot  of  range  as  a function 
of  frequency  for  initial  elevation  angles  of  3°  and  6°  is  shown  in  Figure  1-4. 
For  these  curves  the  frequency  was  sampled  in  0.  25MHz  steps.  The  gaps  in 
the  curves  around  20  MHz  are  due  to  rays  which  are  reflected  once  in  the  F2 
layer  but  penetrate  down  range  on  the  second  oscillation  when  the  F2  layer 
plasma  frequency  is  reduced.  A comparison  of  three  ionospheres  without  an 
ionization  valley  but  with  the  three  different  horizontal  gradients  is  shown 
bel  ow. 

dN 

Case  — — (x=2000  km)  Number  of  Trapped  Rays 


h=2  30  km 

h=185  km 

ID' 

- 0.  0348  el/cc /m 

-0. 0167 

4 

2D' 

-0.  0695 

-0. 0335 

9 

3 d' 

-0. 0978 

-0. 0471 

22 

The  horizontal  gradient  is  given  at  only  one  point  but  its  value  is  indicative 
of  the  relative  behavior  over  a large  region.  For  the  same  vertical  profile 
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the  number  of  trapped  rays  increases  faster  than  linearly  as  a function  of 

the  horizontal  gradient.  No  rays  were  trapped  for  any  of  the  three  cases 

of  the  profile  type  A*  which  differs  from  d'  by  having  a weaker  vertical 

ionization  gradient  in  the  E region  but  a stronger  gradient  in  the  F2  region. 

The  effectiveness  of  ducting  of  rays  in  a two-layer  ionospheres  without  an 

ionization  valley  appears  to  be  influenced  both  by  the  ionization  density  of 

the  underlying  layer  and  the  horizontal  gradients  in  the  region  immediately 

above  the  E layer.  In  TABLE  1-3  the  number  of  rays  ducted  (Nf^  as  a 

function  of  frequency  and  initial  elevation  angle)  is  presented  for  several 

profiles  without  an  ionization  valley  and  for  the  same  horizontal  gradient 

(Case  2).  For  both  sets  of  profiles  the  efficiency  of  long  range  ducting 

increases  with  the  value  of  foE.  Comparing  the  models  that  have  the  same 

value  of  foE  we  observe  that  the  efficiency  is  larger  for  the  profiles  that 

have  both  a region  where  = 0 followed  by  a region  where  (h)  > (h) 

xx  x 

for  h > 160  km,  where  xx  and  x'  denote  the  two  kinds  of  models.  To  deter- 
mine which  condition  produces  the  increase  in  the  ducting  efficiency  we 
traced  rays  through  a model  which  incorporates  the  steeper  vertical  gradient 
without  the  extended  region  of  zero  vertical  gradient.  This  model  is  2D  where 
hmF2=250  km  instead  of  300  km. 

Only  three  rays  were  long-range  ducted  for  this  model  which  is  less 
than  for  both  the  two  previous  models.  Although  a negative  horizontal  gradient 
is  required  to  facilitate  ducting  the  nature  of  its  coupling  with  a vertical 
ionization  gradient  in  a particular  ionospheric  model  is  very  critical.  For 
certain  ionization  profiles  negative  horizontal  gradients  may  represent  a 
necessary  but  not  sufficient  condition  for  efficient  ducting  of  rays. 

For  ionospheres  with  negative  ionization  gradients  the  occurance  of 
valleys  in  the  vertical  profile  appears  to  enhance  the  efficiency  for  ducting. 

In  TABLE  1-4  the  number  of  rays  ducted  is  presented  for  several  ionospheres 
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as  a function  of  both  vertical  and  horizontal  electron  density  gradients. 


TABLE  1-4 

Long  Range  Ionospheric  Ducting 

Number  of  Rays  Ducted 
(Elevation  Angle  s 15°) 


Vertical  Gradients 


Horizontal  Gradients 

I 1 

2 

0 

0 

0 

4 

9 

22 

8 

16 

33 

9 

20 

44 

11 

24 

50 

14 

28 

58 

>62 

>74 

>110 

DCBA  Combination 


For  a given  vertical  profile  model  shows  the  same  behavior  for  the 
different  horizontal  gradients;  i.  e.  Nf^  increases  at  least  as  fast  as  linearly 
with  gradient.  For  a given  horizontal  ionization  gradient  the  ducting 
efficiency  increases  slowly  with  the  deepening  of  the  ionization  valley.  The 
ionospheric  model  DCBA  is  one  in  which  a negative  horizontal  gradient  is 
introduced  along  the  ionization  valley.  The  vertical  ionization  distribution 
has  no  valley  at  the  transmitter  but  one  begins  to  develop  downstream 
reaching  a minimum  at  a range  of  6000  km. 

A ray  pattern  for  model  2DCBA  along  with  the  iso-ionic  contours  is 
illustrated  in  Figure  1-3.  For  two  elevation  angles  of  5"  and  1 (f  , rays  were 
computed  for  frequencies  on  either  side  of  and  including  the  ducting  bands. 
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Examining  the  conditions  encountered  for  EL  = lCf , it  is  seen  that  both  the 
13  MHz-  and  16  Mhz  ray  satisfy  a one-hop  reflection  condition,  returning 
to  the  ground  at  about  D = 1050  and  2400  km  respectively.  Both  the  14 
MHz-ray  and  the  15  MHz-ray  encounter  favorable  injection  conditions  when 
entering  into  the  ionization  valley  and  remain  there,  apparently  tracking 
the  expansion  and  subsequent  contraction  of  the  ionization  valley  with  distance 
with  an  undulation-height  to  undulation- wavelength  ratio  of  about  0.  0275  for 
both  the  14  MHz-ray  (for  example.  Ah  = 32.  88  km  and  X = 1191  km)  and  the 
15.MHz-ray  (for  example.  Ah  = 42.66  km  and  X = 1555  km),  indicating  the 
possibility  that  this  ratio  may  be  invariant  over  the  ducted  frequency  band. 

The  ratio  of  N,  (DCBA)/N.  (A)  varies  between  2.  for  case  3 and  4. 

fa  fa  — 

for  case  but  is  a significantly  larger  change  in  ducting  efficiency  relative 
to  any  of  the  profiles  with  fixed  although  deeper  ionization  valleys.  This 
increase  in  efficiency  holds  even  though  the  gradient  in  the  valley  causes  a 
small  positive  horizontal  ionization  gradient  in  the  first  2000  km  of  the 
path.  An  examination  of  the  ray  trajectories  indicates  this  gradient  may 
aid  in  trapping  some  rays  that  might  otherwise  penetrate  the  E layer  and 
return  to  the  ground  after  one  downward  reflection  from  the  F layer.  This 
suggests  that  from  the  study  of  an  initially  ducted  ray  no  inferences  can  be 
drawn  whether  the  ray  will  or  will  not  continue  to  be  ducted  without  knowledge 
of  the  ionosphere  downstream. 

The  distribution  N^for  the  discrete  sampling  of  the  set  of  ionospheric 
models  of  TABLE  1-3  with  the  largest  horizontal  ionization  gradient  is 
shown  on  Figure  1-5.  Each  successive  ionosphere  shown  includes  all  the 
ducted  rays  of  the  previous  ionosphere  with  a smaller  ionization  valley. 

Since  the  ionosphere  is  stratified  below  110  km  the  boundary  specifying  the 
lower  frequency  limit  of  the  ducting  region  can  be  calculated  analytically 
using  the  secant  law  for  oblique  incidence 
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and 


sin  a 


sin  a 
o 


(*■  + M 

f = foE  l— — -J  sec  (90°  - a) 


(1-5) 


(1-6) 


where  a is  the  elevation  angle  at  the  ground  and  foE  = 3.  5 MHz  and  h =110  km 
o ■k 

for  this  model.  For  the  first  five  ionospheric  models  we  observe  an  increase 
in  the  number  of  rays  trapped  at  any  given  elevation  angle  as  the  ionization 
valley  deepens.  This  appears  to  be  due  to  the  reflection  of  the  ray  on  the 
downward  part  of  the  path  in  the  height  region  between  110  and  160  km  where 
the  ionization  gradient  is  increasing  in  the  direction  of  motion  of  the  ray.  We 
note  that  for  these  ionospheres  there  is  no  horizontal  gradient  in  this  region 
(below  160  km).  For  the  last  case,  3DCBA,  the  vertical  gradient  is  weaker 
in  this  height  region  but  there  is  now  a horizontal  ionization  gradient  which 
enlarges  the  net  positive  gradient  along  the  rays  which  tends  to  reduce  the 
local  elevation  angle  and  return  them  to  the  duct. 


The  number  of  rays  ducted  and  various  horizontal  and  vertical  ionization 
gradients  for  models  with  the  same  foF2  gradients  (Case  2)  is  shown  in 
TABLE  1-5.  Some  of  these  results  are  shown  in  TABLE  1-6  for  each  of  the 
possible  permatations  of  f^,  and  f^.  We  observe  several  examples  (enclosed 
values  in  the  Table)  of  three  layer  monotonic  electron  density  profiles  that 
support  long  range  ducted  modes.  For  the  same  f^  and  foF2  parameters 
increasing  the  underlying  ionization  by  increasing  the  value  of  f results 
in  more  rays  being  ducted. 
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TABLE  1-5 


Long-Range  Ionospheric  Ducting 

Fixed  Horizontal  Electron  Density  Gradient 
(Ionosphere  #2) 


Ionosphere 

Nfq 

dN  . 

2000)  el/cc/m. 

3N 

— (x=2000)  el/cc/m. 

as  1 1° 

ox 

h = 185 

135 

185 

UIl 

135 

100km 

2 BA 

0 

-0. 0076 

0 

4.  27 

0.  88 

3.9 

2CA 

0 

-0. 0076 

0 

4.07 

1. 95 

3.9 

2 DA 

0 

-0. 0076 

0 

3.83 

3.  21 

3.  9 

2A' 

0 

-0. 0335 

-0. 0042 

5.82 

2.  47 

3.  90 

2B' 

0 

-0. 0335 

-0.0042 

5.60 

2.  39 

6.09 

2DB 

2 

-0. 0076 

0 

3.83 

2.  34 

6.  08 

2CB 

2 

-0. 0076 

0 

4.07 

1. 07 

6.  08 

2AA 

2 

-0. 0076 

0 

4.44 

0 

3.  90 

2BB 

3 

-0. 0076 

0 

4.  27 

0 

6.08 

2 DC 

4 

-0. 0076 

0 

3.  83 

1. 27 

3.51 

2C1 

4 

-0. 0335 

-0. 0042 

5.  33 

2.  26 

8.  76 

2AB 

9 

-0. 0076 

0 

4.44 

1 

o 

• 

00 

00 

6.09 

2D' 

9 

-0. 0335 

-0. 0042 

5.02 

2.  13 

11.  93 

2CC 

11 

-0. 0076 

0 

4.07 

0 

3.51 

2BC 

14 

-0.0076 

0 

4.27 

-1.07 

6.  08 

2AC 

15 

-0. 0076 

0 

4.44 

-1. 95 

8.  78 

2D(  DD) 

16 

-0.0076 

0 

3.83 

0 

11.  93 

2C(CD) 

20 

-0. 0076 

0 

4.07 

-1. 27 

11.93 

2B(BD) 

24 

-0. 0076 

0 

4.27 

-2.  34 

11.  93 

2A(AD) 

28 

-0. 0076 

0 

4.44 

-3.  21 

11.  93 

I I 
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TABLE  1-6 


Long-Range  Ionospheric  Ducting 
Number  of  Rays  Ducted 
Fixed  Horizontal  Electron  Density  Gradient 


fF: 

2.  0 

(Ionosphere  #2) 

2.5  3.0 

3.5  MHz. 

i 

2.  0 MHz. 

2 

9 

15 

28 

2.5 

0 

3 

14 

24 

3.  0 

0 

! 2 

11 

20 

3.5 

0 

2 

4 ; 

i 
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The  ionospheres  described  thus  far  have  been  characterized  by  a 
negative  horizontal  ionization  gradient  produced  by  the  value  of  foF2 
decreasing  along  the  direction  of  ray  propagation.  In  Figure  1-6  four  ray 
paths  are  plotted  along  with  the  iso-ionic  contours  for  model  2A*  . For 
the  2?  elevation  angle  two  rays,  separated  by  0.  5 MHz  with  similar  ray 
paths  are  ducted.  In  Figure  1-7  we  show  the  iso-ionic  contours  of  a model 
in  which  foF2  = 8 MHz  over  the  entire  10000  km.  Here  the  negative  hori- 
zontal gradient  is  generated  by  the  change  of  f^  along  the  surface  of  the 
earth.  The  values  of  f^  at  160  km  and  at  spacings  of  2000  km  are  3.  0, 

2.5,  2.0,  2.0,  2.  5 and  3.  0 MHz  and  f is  a constant  3.  5 MHz.  Four  rays 

XL< 

with  the  same  parameters  as  in  Figure  1-6  are  traced  through  this  iono- 
spheric model  and  are  also  shown  in  Figure  1-7.  The  ray  paths  for  the  two 
ducted  frequencies  in  both  models  are  quite  similar.  The  total  number  of 
ducted  rays  for  this  ionosphere  is  40  which  falls  between  the  numbers  for 
models  2 A and  3A  which  are  respectively  28  and  55.  Although  a negative 
horizontal  ionization  gradient  appears  to  be  a necessary  condition  for 
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Ionospheric  Model  with  foF2=constant 


Figure  1-7 


long-range  ducting  the  manner  in  which  this  gradient  is  achieved  is  still 
quite  flexible. 


I.  5 Ducting  Efficiency  and  Analytic  Perturbations  to  Ionospheric  Models 

In  each  of  the  models  considered  thus  far  we  have  calculated  a relative 
efficiency  of  the  model  to  support  long  range  ducted  modes  by  using  a discrete 
sampling  of  communications  circuit  parameters  in  a ray  trace  program.  All 
of  these  models  have  been  characterized  by  a slowly  varying  electron  density 
distribution  whose  horizontal  gradients  have  been  determined  from  profiles 
specified  at  2000  km  intervals.  We  have  also  studied  the  changes  in  the  ducting 
efficiency  of  some  of  these  models  for  various  localized  perturbations  to  the 
electron  density  distributions.  The  nature  of  these  changes  yields  some 
information  on  the  stability  of  long  range  modes  for  a particular  set  of  iono- 
spheric conditions. 

In  order  to  insert  perturbations  to  the  ionosphere  we  have  used  a three 
layer  model  and  vary  the  values  of  f^  for  some  profiles  for  part  of  the  ray 
path.  We  generate  a three  layer  profile  which  is  similar  to  the  two  layer 
model  of  3D'  . Values  of  f^  are  determined  from  the  two  layer  model  at  a 
height  of  160  km  at  each  of  the  profiles  given  along  the  ray  path.  Although 
the  resulting  ionosphere  matches  3d'  at  each  of  the  critical  altitudes  there 
are  deviations  between  the  two  models  at  intermediate  heights  between  110 
and  300  km.  We  now  use  two  sine-squared  profiles  in  this  region  where 
previously  we  only  had  one.  This  results  in  a zero  value  of  the  vertical 
electron  density  gradient  at  160  km.  Both  models  have  non-zero  horizontal 
ionization  gradients  above  110  km. 

Figure  1-8  shows  the  distribution  of  the  48  ducted  rays  for  this  three 
layer  ionosphere  (dotted  area).  For  a given  elevation  angle  these  rays  are 
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Three  Layer  Ionospheric  Model  with  Perturbations 


located  at  frequencies  just  above  those  for  penetration  of  the  layers  at  1 1 0 
and  160  km.  The  22  rays  ducted  for  the  two  layer  3D1  model  were  all 
located  along  the  E layer  frequency- elevation  angle  boundary.  In  the 
three  layer  model  the  number  of  ducted  rays  along  the  E layer  boundary 

* : 

is  reduced  to  5 while  43  rays  are  ducted  at  the  new  FI  region  boundary. 

The  separation  of  the  two  ducting  regions  for  a given  elevation  angle  varies 
between  1.5  and  3.0  MHz.  While  the  negative  horizontal  gradient  enhances 
long-range  ducting  the  prediction  of  circuit  parameters  is  not  only  strongly 
dependent  upon  the  individual  layer  parameters  but  also  on  the  detailed 
nature  of  the  vertical  profile. 

The  ionosphere  that  we  perturb  is  the  three  layer  3D1  model.  The 

disturbance  is  constructed  by  increasing  the  value  of  the  plasma  frequency 

at  h = 160  km  by  y%  at  xq,  y/2%  at  xq±  250  km  and  0%  at  xq±500  km.  In 

this  manner  the  disturbance  is  confined  to  a 1000  km  portion  of  the  path 

and  produces  a change  in  electron  density  centered  at  h = 160  km  that 

vanishes  at  h = 110  and  h = 300  km.  The  results  for  x = 750  km  and  two 

o 

values  of  y,  10%  and  20%  are  also  shown  in  Figure  1-8.  Several  features 
are  apparent  from  the  distributions  for  the  perturbed  and  unperturbed  cases. 

1.  The  amount  of  long  range  ducting  increases.  For  the  10%  per- 
turbation the  number  of  ducted  rays  are  9 and  48  along  the  two  boundaries. 

With  a 20%  perturbation  these  values  are  14  and  44  respectively. 

2.  At  a given  frequency  for  ducting  along  the  FI  layer  the  elevation 
angles  required  are  larger.  For  example  at  20  MHz  the  10%  increase  in 
plasma  frequency  shifts  the  ground  elevation  angle  by  3?  and  the  20%  change 
produces  a shift  of  5°.  This  effect  also  produces  a wider  range  of  available 
elevation  angles  for  some  frequencies.  At  22  MHz  the  interval  of 
elevation  angles  changes  from  3°  to  5.5°  to  8°  for  the  three  ionospheres. 
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3.  The  maximum  frequency  for  ducting  is  reduced  by  0.5  MHz. 

This  analysis  was  also  carried  out  for  the  same  values  of  y but  with 

x = 500  km  (250  km  closer  to  the  transmitter).  The  parameters  for  the 
o 

ducted  rays  along  the  E layer  were  shifted  but  the  number  of  them  did  not 
increase.  This  is  probably  due  to  the  discrete  sampling  of  a very  narrow 
region  of  parameters  along  this  boundary.  The  number  of  ducted  rays  along 
the  FI  boundary  increased  to  46  for  both  perturbations.  These  showed  the 
same  tendency  of  narrowing  the  available  frequency  interval  and  widening 
the  angular  interval  as  the  previous  case  of  xq  = 750  km. 

We  conclude  that  while  relatively  small  changes  from  a given 
plasma  frequency  distribution  do  not  destroy  the  possibility  of  supporting 
long-range  modes  they  can  substantially  change  the  required  circuit  param- 
eters. 

Our  sampling  interval  revealed  no  ducted  modes  in  either  a two  or 

three  layer  stratified  ionosphere.  Ducting  was  observed  when  negative 

horizontal  gradients  were  applied  over  a large  portion  of  the  total  ionospheric 

cross  section.  It  is  of  interest  to  know  whether  these  gradients  need  exist 

over  such  a large  portion  of  the  potential  path  or  are  only  required  near  the 

transmitter  in  order  to  obtain  injection  into  a ducted  mode.  We  have  used 

the  method  described  above  to  introduce  gradients  over  a restricted  region 

of  a stratified  ionosphere.  We  have  considered  two  basic  models  with  fixed 

parameters,  f = 3.  5,  f =5.  0,  f = 8.  0 MHz  and  f = 2. 5,  f = 5. 0,  f,  = 8.  0 

MHz.  Both  of  these  models  were  subjected  to  perturbations  with  y = 20% 

and  x = 500  km.  No  ducted  rays  were  observed  for  the  model  with  f = 2.5 
o E 

MHz.  Three  rays  along  the  E region  frequency- elevation  angle  boundary 
were  ducted  for  the  model  with  f = 3.  5 MHz.  The  elevation  angles  for 
these  rays  were  between  4.5°  and  7.(f  which  corresponds  to  the  rays 
intersecting  the  bottom  of  the  ionosphere  (h  = 90  km)  at  ground  ranges  of 
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640-540  km.  The  presence  of  a negative  horizontal  gradient  in  this  region 
is  sufficient  to  inject  into  a duct  even  though  the  ionosphere  remains 
stratified  for  the  remaining  9000  km  of  the  path. 


I.  6 Ducting  Efficiency  and  Random  Perturbations  to  Ionospheric  Models 

The  results  of  the  previous  section  are  helpful  in  accessing  the  changes 
in  ducting  efficient  due  to  differences  in  actual  ionospheric  parameters  from 
those  predicted  by  a median  model.  These  may  represent  measureable  but 
unknown  differences  over  a large  region  of  space,  e.g.  200-2000  km.  It  is 
also  of  interest  to  make  estimates  about  the  stability  of  various  ducted  modes 
to  more  random  and  highly  localized  fluctuations  of  ionospheric  conditions. 

We  have  developed  a scheme  which  allows  one  to  combine  the  modeled  effects 
of  random  scatterings  with  a deterministic  ray  trajectory.  Our  two-dimensional 
ray  tracing  program  calculates  the  ray  coordinates  at  a discrete  series  of  points 
along  the  surface  of  the  earth.  We  solve  for  h^(x^)  and  a^(x^)  for  equal  spacing  in 
x..  In  our  scheme  the  ray  is  repeatedly  subjected  to  a series  of  perturbation 
using  a scale  size  that  is  related  to  the  probability  of  the  ray  being  disturbed. 
When  the  ray,  advancing  in  10  km  steps  reaches  a distance  corresponding  to 
the  average  local  scale  size  the  ray  direction  is  perturbed. 

The  magnitude  and  sign  of  the  angular  perturbation  are  derived  from  a 
random  Gaussian  distribution.  Both  the  scale  size  and  the  angular 
perturbation  Acn  are  weakly  dependent  upon  the  electron  density  through  the 
following  models: 


where  N is  the  electron  density  per  cubic  centimeter,  Nq  = 10  el/cc  and 

both  a and  b have  been  chosen  to  be  0.  1 in  this  study.  S and  a are  chosen 

o o 

to  describe  the  nature  of  the  fluctuations  present  in  a particular  ionosphere. 

We  have  used  various  combinations  of  S , or  with  100  km  sS  5 400  km  and 

o o o 

0.005  radians  £ a £ 0.02  radians.  In  the  expression  for  the  angular  pertur- 
o 

bation  g is  a Gaussian  random  number  with  average  zero,  a = 1.0  and  g is 
limited  to  values  in  the  interval  to  ± 6ff.  Since  the  value  of  changes  with 
each  step  in  the  integration  we  use  an  averaged  value  of  the  scale  length 
over  several  steps  in  order  to  determine  when  to  apply  the  angular  pertur- 
bation. When  this  average  value  is  exceeded  the  ray  trace  continues  with 
the  value  of  a.  replaced  by  a.  + La^.  The  average  angular  change  in  a 10  km 
step  in  the  deterministic  ray  trace  calculation  is  approximately  0.  7? . Depending 
upon  the  perturbation  parameters  chosen  the  instantenous  changes  due  to  the 
scattering  are  less  than  4°. 

We  have  used  the  model  3D*  to  determine  the  effect  of  this  repeated 
scattering  process.  This  was  first  carried  out  for  an  unducted  ray  with 
f = 21  MHz  and  a = 10°  launched  from  the  ground.  The  unperturbed  ray 
penetrates  the  E layer  maximum  height,  is  reflected  in  the  F2  layer  and 
has  a ground  range  of  2600  km.  The  perturbation  scheme  was  applied  100 
times  using  the  same  initial  conditions  for  a given  ray. 

The  analysis  was  performed  for  three  sets  of  values  of  Sq  and  aQ.  The 
frequency  distribution  for  the  calculated  ranges  for  each  of  these  cases  is 
shown  in  Figure  1-9.  The  second  and  third  distributions  each  have  one 
parameter  in  comnion  with  the  first  case  which  has  S0  = 200  km  and 
aQ  = 0.01  radians  so  that  we  can  isolate  the  effect  of  changing  each  parameter 
on  the  resultant  distribution.  For  the  second  distribution  the  scale  length  is 
unchanged  but  we  one-half  the  value  of  the  average  scattering  angle  while  for 
the  third  case  the  scattering  angle  remains  unchanged  but  we  double  the  scale 
length  (decrease  the  number  of  applied  perturbations  by  a factor  of  two).  We 
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Random  Perturbation  of  Ray  Direction 


Range  (km.  ) 

frequency  = 21.0  MHz  Elevation  angle  = 10. 


observe  that  for  this  algorithm  the  magnitude  of  each  scattering  has  a considerably 
stronger  effect  on  the  range  distribution  than  the  number  of  such  perturbations. 

We  have  sampled  the  3D'  for  an  elevation  angle  of  3°  every  0.  1 MHz  in 
the  vicinity  of  the  penetration  of  the  E layer  maximum.  We  found  that  for 
f = 18.  3 MHz  the  ground  range  was  2000  km  and  the  long-range  ducted 
region  extends  from  18.4  - 19.0  MHz.  We  have  chosen  to  study  the  effect 
of  repeated  small  angle  scatterings  on  three  rays;  18.  3 MHz  (an  unducted 
ray),  18.  5 MHz  (a  ducted  ray  very  close  to  the  edge  of  the  frequency  band 
for  ducting),  and  18.  7 MHz  (a  ducted  ray  in  the  center  of  the  ducted  band). 

Ten  perturbed  rays  for  f = 18.5  MHz,  experiencing  approximately  ten 
angular  perturbations  per  2000  km,  are  illustrated  in  Figure  I- 10.  For  the 
chosen  perturbation  parameters  we  observe  that  some  rays  remain  in  the 
elevated  duct  while  others  are  deflected  to  the  ground.  As  a measure  of  the 
ducting  stability  we  use  the  percentage  of  the  total  number  of  rays  whose 
range  remains  greater  than  5000  km.  The  results  for  the  three  ray 
frequencies  chosen  and  various  combinations  of  perturbation  parameters 
are  given  in  TABLE  1-7. 

Keeping  in  mind  that  the  values  in  the  table  are  effected  by  the  total 
number  of  rays  sampled  for  each  set  of  parameters  some  consistent 
observations  can  be  made.  For  the  frequency  that  is  not  ducted  in  the 
unperturbed  ionosphere  the  introduction  of  the  scattering  mechanism 
produces  some  long-range  rays.  The  probability  of  entering  the  long-range 
duct  is  seen  to  increase  with  the  frequency  and  size  of  the  scatterings. 
Correspondingly  the  initial  conditions  which  normally  yield  long-range  rays 
are  observed  to  suffer  some  losses  due  to  the  repeated  scatterings.  For 
both  rays  the  probability  of  leaving  the  duct  also  increases  with  the 
frequency  and  size  of  the  scatterings.  The  stability  for  ducting  does  not 
appear  to  be  strongly  correlated  with  the  position  of  the  frequency  in  the 
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ducting  band.  For  the  last  three  sets  of  scattering  parameters  the  stability 
of  ducting  is  larger  for  the  frequency  which  is  in  the  middle  of  the  ducting 
band  than  that  for  the  frequency  that  is  at  the  edge  of  the  band  while  for  first 
two  sets  of  parameters  the  results  are  reversed. 


I.  7 Elevated  Sources  and  Injection  and  Maintenance  of  Ray  Ducts 

We  have  had  several  indications  that  the  conditions  necessary  for  the 
support  of  long-range  ducted  modes  from  the  ground  may  be  separated  into 
two  parts;  those  required  for  injection  into  an  elevated  ionospheric  duct  and 
those  required  to  maintain  the  ray  in  a duct  once  it  has  been  injected.  A small 
perturbation  to  a stratified  ionosphere  in  a localized  region  was  sufficient 
to  provide  the  means  for  injection  when  the  underlying  ionization  was  above 
some  value.  Once  injected,  this  stratified  monotonically  increasing  profile 
was  sufficient  to  maintain  the  duct  for  9000  km.  For  a series  of  models  all 
with  the  same  horizontal  gradient,  the  introduction  and  deepening  of  a valley 
in  the  electron  density  profile  increased  the  efficiency  of  the  ducting. 

In  order  to  further  examine  these  two  conditions  we  have  traced  rays 

through  ionospheric  models  with  elevated  sources.  Figure  I- 11  shows  the 

distribution  N,  for  long-range  rays  for  model  3A  with  the  source  at  160  km. 
lot 

The  search  is  at  4°  intervals  in  local  elevation  angle  and  every  2 MHz  in 
frequency.  Figure  1-12  is  a plot  of  ray  paths  for  a = 0°  and  f between  12  and 
34  MHz.  We  observe  many  more  ducted  modes  for  a source  placed  at  the 
ionization  minimum  than  for  one  at  the  ground.  For  one  of  these  modes  to 
be  reached  from  the  ground  the  underlying  ionization  and  the  horizontal 
gradients  have  to  be  such  that  if  the  ray  reaches  160  km  it  will  have  local 
elevation  angle  which  falls  within  the  boundary  of  the  curve  shown  in  Figure 
1-11. 
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The  changes  in  the  N,  distribution  for  an  elevated  source  brought  about 

fa 

by  a horizontal  ionization  gradient  is  shown  in  Figure  1-13  where  the  solid 
curve  is  the  boundary  for  the  model  3A  and  the  dashed  curve  is  for  a 
stratified  model  with  f = 3.5  MHz,  f,  = 2.0  MHz  and  f_  =8.0  MHz.  This 
latter  model  does  not  support  any  ducted  modes  for  a ground  based  trans- 
mitter. Only  positive  angles  are  traced  and  the  region  of  interest  is  shown 
finer  detail  with  a 0.  5 MHz  and  0.5°  sampling  interval.  At  small  elevation 
angles  the  effect  of  the  horizontal  gradient  is  to  reduce  the  maximum 
frequency  that  will  support  a ducted  modes.  At  elevation  angles  above  lCf 
the  horizontal  gradient  tends  to  increase  the  maximum  possible  ducted 
frequency.  Also  shown  in  Figure  1-13  are  the  curves  representing  the 
transformation  between  ground  elevation  angles  and  frequencies  to  elevation 
angles  at  h = 160  km  and  the  same  frequency  obtained  by  tracing  rays  up  to 
a height  of  160  km.  This  height  is  reached  at  various  points  downstream 
from  the  transmitter.  They  can  not  be  made  to  correspond  with  the  boundary 
curves  for  the  elevated  source  since  they  all  originate  at  the  origin  of  the 
coordinate  system.  If  an  exact  matching  of  elevation  angle  and  the  corre- 
sponding ground  range  could  be  made  the  overlap  of  the  two  curves  would 
denote  all  the  long-range  ducted  modes  that  are  possible.  The  cross  hatched 
region  shows  which  rays  are  in  fact  ducted  for  the  3A  ionospheric  model. 

The  horizontal  gradient  enhances  the  overlap  area  between  the  two  regions. 
While  this  method  shows  the  nature  of  the  required  gradients  for  injection 
it  is  computionally  impractical  to  use  for  different  ionospheric  models.  The 
case  considered  here  used  a model  with  no  horizontal  gradients  below 
h = 160  km  to  calculate  the  required  transformation  of  elevation  angles. 


II.  THREE  DIMENSIONAL  ELECTRON  DENSITY  MODE LS 


II.  1 Introduction  in  the  Polar  Region 

In  the  previous  section  we  have  used  simple  electron  density  models 
in  the  study  of  a particular  HF  communication  problem.  In  order  to  study 
this  problem  further  and  for  other  ionospheric  applications  we  have  con- 
tinued our  work  in  the  development  of  three  dimensional  electron  density 
models  in  the  polar  region.  In  reports  ^ ^ ^ we  have  described  the 
two  elements  of  such  models.  They  consir*  of  empirical  models  for 
individual  ionospheric  layer  parameters  and  a method  for  generating  a 
continuous  vertical  electron  density  distribution  from  these  parameters. 
These  vertical  profiles  have  involved  the  use  of  simple  analytic  functions 
which  maintain  continuity  of  both  the  electron  density  and  its  height  de- 
rivative. Attempts  have  been  made  to  use  these  one  dimensional  dis- 
tributions to  construct  physically  accurate  three  dimensional  ionospheric 
models.  The  extension  in  spaci  is  made  using  numerical  interpolation 
techniques  between  the  predicted  vertical  profiles  at  specified  points  on  a 

Examination  of  the  resulting  isoionic  contours  in  two  dimensions  of 
this  model  (POLAR  76)  revealed  several  problems  and  deficiencies  with 
the  methods  used  in  constructing  the  model.  We  have  constructed  a new 
vertical  electron  density  model  (POLAR  77)  in  an  attempt  to  improve  the 
model.  No  changes  were  made  to  the  interpolation  procedures  used  to 
extend  the  profiles  in  space.  A second  order  interpolation  method  is  used 
and  appears  to  be  adequate  to  describe  all  possible  horizontal  electron 
density  gradients. 

In  the  development  of  POLAR  77  we  have  also  found  it  necessary  to 
modify  some  of  the  models  for  the  individual  ionospheric  parameters  so 
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that  the  electron  density  distributions  presented  here  using  the  POLAR  76 
model  are  not  the  same  as  would  be  constructed  from  the  model  of  earlier 
report.  We  have  also  modified  the  daytime  profile  of  this  model  such  that 
it  always  is  made  up  of  cubic  equations  between  the  E and  FI  layers  and  the 
FI  and  F2  layers  independent  of  the  relative  values  of  the  individual  layer 
parameters.  This  removes  some  unphysical  gradients  that  were  predicted 
when  the  zenith  angle  is  in  the  region  75°  s y s 90^ . Both  models  use  the 
same  parameters  based  upon  median  ionospheric  conditions  when  the  following 
are  specified:  latitude,  longitude,  universal  time,  sunspot  number  and  index. 

II.  2 Development  of  Electron  Density  Model 

The  two  electron  density  models  are  identical  in  the  height  region 

above  h F2  and  in  the  daytime  for  altitudes  below  h E.  The  models 
m m 

differ  in  the  region  between  h E and  h F2.  The  problems  that  arose 

m m 

with  POLAR  76  in  this  region  include  the  following: 

a)  The  disappearance  of  individual  layers  as  a function  of  space 
and  time  and  the  resulting  discontinuity  due  to  the  use  of 
different  analytic  functions  depending  upon  the  number  of 
layers  present,  (e.g.,  the  solar  E and  FI  layers). 

b)  The  difficulty  in  predicting  the  extrema  position  and  magnitude 
of  the  electron  density  in  the  region  between  100  and  200  km 
when  only  the  auroral  E and  F2  layers  are  present  at  widely 
separated  heights. 

c)  Changes  in  the  shape  of  the  vertical  profile  when  only  the 
auroral  E and  F2  layers  are  present  and  the  sign  of  the 
difference  between  the  magnitude  of  foF2  and  foE  changes  at 
locations  which  are  close  to  each  other  in  space  and  time. 
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In  order  to  eliminate  some  of  these  problems  we  have  replaced  the  various 
analytic  functions  used  in  the  vertical  profile  model  with  a single  set  of 
functions  containing  all  the  possible  ionospheric  parameters.  The 
coefficients  of  terms  representing  individual  layers  were  constructed  in  a 
manner  which  attempts  to  maintain  continuity  in  the  other  spatial  variables. 
This  required  the  extension  of  the  individual  models  of  various  parameters 
into  regions  where  they  did  not  exist.  An  example  of  this  is  the  prediction 
of  an  electron  density  at  a height  corresponding  to  hmFl  for  y (zenith  angle) 

> 9(f. 

In  order  to  describe  the  electron  density  in  the  region  between  hmE 
and  hmF2  we  make  use  of  median  predictions  for  the  following  ionospheric 
parameters: 

foF2,  hmF2,  ymF2,  foFl,  hmFl,  ^0^s0^ar 

hmE  , foE  , and  hmE  , 

solar,  auroral,  auroral 

We  first  assumed  that  the  electron  density  could  be  described  by  the  following 
fifth  order  polynomial: 

5 4 3 2 

(1)  N (x)  = ax  + bx  + cx  + dx  + ex  + f where 
m 

x = h-  hmE 


hmF2-  hmE 


In  order  to  determine  the  six  coefficients,  we  used  the  above  set  of  iono- 
spheric parameters  to  specify  the  electron  density  and  its  height  derivative 
at  three  points  for  OSxSl.  Using  this  model,  all  the  coefficients  were 
constrained  to  be  continuous  functions  of  latitude,  longitude  and  universal 


41 


f 


time.  This  ensured  that  no  gradients  existed  in  the  three-dimensional 

model  that  were  not  consistent  with  those  that  existed  for  the  median  models 

of  the  ionospheric  parameters.  Since  the  diurnal  variation  in  some  of  the 

3 

parameters  (e.g.  foE)  varied  by  a factor  of  10  , we  also  used  the  expression 
log  N^(x)  in  the  polynomial.  This  form  produced  more  realistic  height 
gradients  in  the  profiles. 

In  general,  the  polynomial  (1)  has  four  extrema  and  at  least  one  real 

root  (i.  e.  N (x)  = 0 for  some  value  of  x).  Three  of  the  extrema  are 
m 

specified  by  the  parameters.  In  order  for  the  resultant  profile  to  be 
physically  valid,  the  fourth  extrema  and  any  real  roots  cannot  occur  in  the 
interval  OfixSl.  This  proved  to  be  a problem  as  the  polynomial  obtained 
from  the  three  data  often  predicted  a spurous  maximum  in  the  height  region 
of  interest. 

In  order  to  remove  this  difficulty,  we  used  a fourth  order  polynomial 
of  the  form: 

4 3 2 

(2)  log  N (x)  = ax  + bx  + cx  + dx  + e 


The  ionospheric  parameters  were  used  to  form  the  following  constraints: 


(3) 

(4) 

(5) 

(6) 

(7) 


d_ 

dx 


d_ 

dx 


N (0)  = 
m 


log  Nm(0) 


N (1)  = 
m 


log  Nm(l) 


N E 
m 

= 0 

N F2 
m 

= 0 


N (x. ) = N,  where  x,  = h,  -hmE 
m 11  1 _1 

hmF2-  hmE 


42 


In  order  to  reduce  the  computational  complexity  we  introduce  the  cubic 
equation  which  satisfies  conditions  (3)  - (6) 


P(x)  = -2  (l°g  NmF2  * loS  NmE)  x3  + 3 (log  NmFZ  - log  NJ)  x2 


m 


+ log  N E 
m 


The  resulting  fourth  order  polynomial  is 


"&>*  w - pM' 


l»g  n„w  = P(x)  + { — 5 — ‘"H  (*2(x-  l)2) 


2 / IS2 

Xj  (Xj-1) 


II.  3 Layer  Parameters 

For  solar  daytime  conditions  all  the  parameters  that  are  necessary  to 
define  the  polynomial  are  available  from  the  individual  models.  We  have 
constructed  models  for  nighttime  conditions  and  modified  the  definition  of 
the  parameters  given  in  (3)  - (7)  so  that  they  are  continuous  in  space  and 
time. 


E Layer 

The  solar  E layer  parameters  foE  and  hmE  are  described  in  reference 

s s 

4.  We  have  modified  the  definition  of  hmE 

s 


hmE 

s 

= 100. 

+ 20 

hmE 

s 

= 100. 

+ 10 

loge(secy)  x<87.1° 

log  (sec  87.1°)  = 160.  km  y ^ 87.  1° 
e 
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For  X > 9C f 


3 

N^Eg  = 2x10  + 30.  0 S el/cc  where  S is  the  average 

sunspot  number  and  h E =110.  km.  The  southern  and  northern 

m s 

boundaries  of  the  auroral  oval  are  defined  by 

Magnetic  Day  = 78.  0 - 2.  5 deg.  6.0st<sl8.0 

Magnetic  Night  $ = 71.0  - 2.  5 K cos  (t  - 1)),  t <6.0  or  t > 18.0 

s p \1Z  cic  c 

and  + Kp  + 4. 0 where  t is  the  corrected  geomagnetic 

time  and  X is  the  corrected  geomagnetic  latitude.  When  $ £ X s t , the 
m ot,  s m N 

auroral  oval  blanketing  frequency  is  calculated  using 

fA  = 1. 3 + 0.  3 (4.  0 + 0. 05  K ) 

A p 

For  solar  night  this  is  replaced  by 


+ 


,2)1/2 


The  final  E layer  parameters  foE  and  hmE  are  generated  according  to  the 
rules 


(foE)4  = (foE  )4  + |(foEA)4 

S tL  A 

hmE  = (foE  )2  hmE  + (foE  )2  hmE 

S S A A 

(foEs)2  + (foEA)2 


F2  Layer 

The  calculation  of  the  maximum  density  N F2  and  the  semithickness 

m 

of  the  F2  layer  are  given  in  reference  4 and  5.  The  model  for  hmF2  has 
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been  slightly  modified.  The  seasonal  variation  of  h^F2  is  incorporated 


in  the  behavior  of  noon  and  midnight  values. 

2 

h F 2 (noon)  = 197.0  + 0.  79S  - 0. 0011S 
m 

+ 0. 55(X  - 45.  CP) 

m 

h F2  (midnight)  = 297.0  + 0.603S  + 0.55(X  - 45.  0°) 

m m 

The  diurnal  variation  of  h F2  is  divided  into  three  separate  cases.  In  the 

m 

first  case  there  is  no  diurnal  variation.  This  occurs  when  the  sun  does  not 
set  or  when  the  sun  does  not  rise. 

h^F2  = h^F2  (noon)  if  < 1 02f* 

h F2  = h F2  (midnight)  if  y(noon)>102f 
m m 

where  \ is  the  local  zenith  angle.  The  remaining  two  cases  depend  upon  the 
value  of  x(noon)  and  the  behavior  of  h^F2  is  shown  in  the  two  following 
curves.  Linear  interpolations  for  hmF2  are  performed  from  one  of  the 
two  graphs  for  a given  value  of  In  order  to  determine  which  case 
applies,  the  four  parameters;  x(noon),  x(midnight)»  T(x  = 87°)  are  calculated 
for  the  given  day  of  the  year. 
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FI  Layer 

A twenty-four  hour  model  for  an  electron  density  at  a height  between 
the  E and  F2  layers  is  constructed  using  the  solar  FI  layer  and  the  various 
parameters  of  the  F2  layer.  The  model  for  N^F^hj)  for  y<90°  is  de- 
scribed in  reference  1.  The  height  of  this  intermediate  layer  corresponds 

to  h FI  and  is  given  by 
m 

hj  = hmFl  = 156.  + 0.15S  + 45.  log^sec  x),y  £ 7<f 

h,  = 156  + 0.  15  S + 45.  log  (sec  7(f) 

1 e 

= 204.  3 + 0. 15  S X>7(f 

For  y>9  (f 

N,  = N FI  = 103  + 70.  OS  el/cc. 

1 m 

A cubic  equation  is  used  with  N F2,  hmF2  and  ay  F2,  where  a is  a 

m m 

parameter  between  1.25  - 2.  25,  to  calculate  the  electron  density  at  h.  that 
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would  be  predicted  by  an  F2  layer  alone.  This  value,  N^F21  is  combined 

with  N FI  to  produce  an  anchor  point  on  the  profile  at  all  times  of  the  day. 
m 

The  weighting  procedure  is  the  following: 
h F2  - h 

A,  s —22 — with  A,  always  <1.0, 

1 a y F2  1 

m 


2 2 
A2  - 1.0-  A, 

struction  of  the  profile  is  given  by 


and  the  density  at  h^  used  in  the  con- 


N,  = N FI  A?  + N F21  A^. 
1 ml  m l 


As  h,  becomes  more  removed  from  h F2  the  less  influence  the  F2  layer 
1 m 

has  on  the  determination  of  N^.  At  night  A^  -*  1 and  approaches  the  value 

N FI  (x>9Cf).  During  the  day  for  periods  close  to  the  equinox  the  weighting 
m 

factors  are  of  the  same  size  and  the  two  densities  are  close  to  the  same 
value.  The  procedure  then  produces  a smoothed  value  for  for  two  inde- 
pendent parameter  models.  For  periods  close  to  the  solstice  h^F2  and  h^ 
can  be  approximately  the  same  (as  predicted  from  the  models).  In  this 
case  the  F2  layer  becomes  the  important  contributor  to  N^.  This  procedure 
reflects  our  assessment  of  the  relative  merits  of  the  individual  empirical 
models  for  the  FI  and  F2  layer  parameters. 


II.  4 One  and  Two  Dimensional  Electron  Density  Distributions 


Table  II- 1 contains  the  ionospheric  parameters  for  various  layers  for 
several  locations  along  a great  circle  path  between  California  and  England 
for  April  15  at  UT  = 2100  hours.  Two  of  the  sets,  0 = ( f and  75°,  are  during 
solar  day  and  night  respectively  while  the  remaining  three  locations  are  close 
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to  the  terminator  at  sunset.  In  figure  II- 1 we  show  the  dependence  of  the 

POLAR  77  model  using  these  parameters  for  various  values  of  a.  Three 

profiles  are  shown  with  a = 1 . 25 , 1.75,  and  2.  25  along  with  the  profile 

constructed  using  the  POLAR  76  model.  It  is  not  our  intent  to  compare 

one  model  with  another  but  to  determine  under  what  conditions  either  can 

be  used  to  describe  realistic  median  ionospheric  conditions.  For  daytime 

conditions  both  POLAR  76  and  POLAR  77  with  a = 1. 25  and  1. 75  contain 

the  general  features  of  the  three  layers.  The  height  derivative  for 

on 

POLAR  76  is  more  slowly  varying  in  a larger  region  around  the  layer 
heights.  POLAR  77  with  a = 2.  25  produces  an  unphysical  peak  in  the 
electron  density  larger  than  N^F2  at  a height  between  h^  and  hmF2. 

If  we  restrict  our  examination  of  the  nighttime  profile  to  > 7000  el/cc 
(0.  75  MHz.  ) we  observe  that  the  gradient  decreases  with  increasing  values 
of  c.  For  a = 2.  25  the  profile  predicts  relatively  large  densities  in  the 
height  region  between  160  and  200  km.  The  POLAR  76  model  is  a single 
parabolic  layer  with  a semithickness  of  1.414  y^F2.  For  a fixed  value  of 
a (1.75)  we  wish  to  examine  the  relative  spatial  and  temporal  dependence 
of  the  two  models. 

Figure  II- 2 shows  the  profiles  using  POLAR  77  for  the  three  locations 
bracketing  the  terminator  and  one  for  low  zenith  angle.  Figure  II- 3 also 
contains  the  same  profiles  using  POLAR  76.  In  the  height  region  between 

h.  and  h F2  the  two  models  are  almost  identical  near  midday  and 

1 m dN 

remain  similar  until  past  sunset  when  -rr—  becomes  larger  with  the  POLAR  77 

oh 

model.  In  the  height  region  below  h^  the  POLAR  77  model  predicts  a 
smaller  value  of  the  total  electron  content  (TEC)  for  low  zenith  angle  and 
a larger  TEC  during  solar  night.  The  shape  of  the  profiles  for  the  POLAR  77 
model  changes  more  slowly  than  those  of  the  POLAR  76  model  as  the 
terminator  is  crossed. 
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Figures  II-3  through  LL-6  are  iso-ionic  contour  maps  of  electron 

density  for  a 75°  great  circle  path  with  the  origin  at  45° N,  114°W  and  an 

azimuthal  angle  of  38.2°  East  for  the  time  and  date  given  in  TABLE  11-1. 

The  height  range  is  between  80.  and  400.  km  and  the  contour  labels  are 

_5 

electron  density  x 10  . The  origin  is  during  solar  day  and  the  path  is 

such  that  we  cross  the  auroral  oval  at  approximately  25°  and  exit  the  oval 
and  cross  the  terminator  between  45°  and  5tf  . The  contours  are  con- 
structed using  a quadratic  smoothing  technique  from  profiles  that  are 
given  every  3f*  in  arc  length  and  every  10.0  km  in  height.  Figure  II-7  is 
plot  showing  the  boundaries  of  the  auroral  oval  and  the  terminator  for  this 
path  as  a function  of  universal  time.  These  boundaries  are  contained  in 
the  individual  models  of  the  ionospheric  parameters.  Figure  11-3  is  con- 
structed using  the  POLAR  76  model  and  figures  11-4  through  II- 6 using  the 
POLAR  77  model  with  a - 1.  25,  1.75,  and  2.25.  For  heights  above  h^F2 
all  the  models  are  identical.  Several  features  of  the  nighttime  part  of  the 
ionosphere  are  clearly  apparent.  The  POLAR  76  model  and  the  POLAR  77 
model  with  a = 1. 75  are  very  similar.  The  POLAR  77  model  with  a = 2.  25 
produces  significantly  larger  electron  densities  at  lower  altitudes  than  the 
others  and  the  model  with  a = 1. 25  produces  an  exceptionally  thin  F2  layer. 
This  leads  to  very  large  horizontal  gradients  at  the  terminator.  Both 
POLAR  76  and  POLAR  77  with  a = 2.  25  have  a large  horizontal  gradient 
in  150-200  km  region  at  the  auroral  oval  boundary.  The  other  two 
POLAR  77  models  generate  a smoother  transition  between  the  auroral 
E layer  and  the  corresponding  F layers.  Figures  II-8,  POLAR  76  and 
II-9,  POLAR  77,  a = 1.75  are  the  same  as  figures  II-3  and  II-5  but  the 
profiles  are  generated  at  1°  intervals  in  arc  length.  In  the  region  of  the 
terminator  the  horizontal  electron  density  gradient  does  not  change 
significantly  with  the  POLAR  76  model  in  the  height  region  just  below 
hmF2.  The  corresponding  gradient  in  the  POLAR  77  model  shows  a 
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Contours  of  Constant  Electron  Density 


GREAT  CIRCLE  DISTANCE 


Contours  of  Constant  Electron  Density  (x  10  el/cc) 


GREAT  CIRCLE  DISTANCE 
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sharper  increase  in  this  region. 

Figures  II- 10,  POLAR  76  andll-ll,  POLAR  77,  a=  1.75  are  iso-ionic 
contours  constructed  at  3°  intervals  for  the  same  path  and  date  but  for 
UT  = 400  hours.  As  can  be  seen  from  Figure  II- 7 this  path  is  entirely 
during  the  solar  night  but  crosses  the  auroral  oval  twice.  In  these  two 
spatial  regions  the  profiles  are  no  longer  monotonically  increasing  up  to 
hmF2.  In  the  POLAR  76  model  the  oval  has  the  additional  effect  of  pro- 
ducing a larger  electron  density  at  heights  below  the  F2  peak  whereas  in 
the  POLAR  77  model  the  effect  of  the  oval  is  almost  entirely  contained  in 

the  region  close  to  h E. 

m 

The  development  of  the  POLAR  77  model  has  revealed  deficiencies  in 
the  individual  ionospheric  parameter  models.  The  changes  made  to  these 
parameter  models  have  resulted  in  an  improvement  of  both  the  profile 
models.  Both  models  have  proved  to  be  useful  in  three  dimensional  ray 
tracing  applications.  Further  improvement  in  predicted  vertical  electron 
density  distributions  would  be  based  upon  more  accurate  ionospheric 
parameter  models.  The  modeling  of  the  F2  and  FI  layers  from  observa- 
tions at  sunrise  and  sunset  cannot  be  performed  separately.  The  ionization 
in  the  region  between  200-250  km  cannot  be  separated  into  individual  layers  since 
the  resulting  FI  and  F2  models  will  not  be  internally  consistent  with  each  other. 

When  hj  = hmF2  the  POLAR  77  model  does  not  have  sufficient  data  to 
produce  a profile  in  the  region  between  h^E  and  h^F2.  Under  these 
conditions  the  POLAR  76  model  would  predict  a zero-width  F2  layer. 

Problems  such  as  these  reveal  the  inadequacy  of  constructing  profiles 
based  upon  specifically  identifiable  properties  of  individual  ionospheric 
layers.  When  the  description  of  the  electron  density  in  terms  of  layers 
breaks  down  the  profile  model  also  breaks  down. 
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Contours  of  Constant  Electron 
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Contours  of  Constant  Electron  Density 
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III.  EARTH  DETACHED  RF  PROPAGATION 
INDUCED  BY  FIELD  ALIGNED  IRREGULARITIES 


There  are  modes  of  radio  propagation  controlled  by  the  regular  iono- 
spheric regions  by  which  HF  and  even  VHF  radio  waves  can  travel  to  great 
distances  with  low  attenuation  without  intermediate  ground  reflections.  An 
important  class  of  these  modes  is  the  Earth  Detached,  or  Chordal  Hop, 
mode  whereby  a ray  is  refracted  in  the  ionosphere  so  as  to  travel  above 
the  earth's  surface,  reentering  the  ionosphere.  Such  modes  of  propagation 
could  be  extremely  useful  for  long  range  HF  communication  or  over  the 
horizon  backs catter  radar  applications. 

However,  it  is  clear  that  in  a spherically  symmetric  ionosphere  such 
modes  cannot  be  excited  by  a ground  based  transmitter.  On  the  other  hand, 
the  naturally  occurring  horizontal  gradients  in  the  ionosphere  could  create 
sufficient  ionospheric  "tilt"  to  inject  a ray  into  earth  detached  propagation. 
Generally,  even  at  grazing  angles,  a ray  leaving  the  ground  arrives  at  iono- 
spheric altitudes  with  too  large  an  angle  with  the  layer  for  injection  into 
earth  detached  propagation.  Furthermore,  the  occurrence  of  favorable 
conditions  for  such  injection  is  so  fortuitous  and  for  such  short  duration  as 
to  preclude  reliance  on  the  mode  for  routine  use. 

Other  natural  mechanisms  capable  of  coupling  HF  radio  waves  into  or 
out  of  earth  detached  propagation  consist  of  (a)  sporadic  E ionization,  (b) 
auroral  columns  and  curtains,  (c)  meteor  trails,  (d)  equatorial  F- region 
ionization  troughs,  and  (f)  travelling  ionospheric  disturbances.  Clearly, 
these  all  lack  the  reliability  criterior  for  routine  utilization  of  the  mode. 

Artificial  scattering  centers  such  as  rocket- launched  ionization  trails, 
rocket-launched  ionospheric  seeding,  or  an  RF  induced  ionization  cloud, 
are  also  conceivable  and  have  been  shown  to  be  effective  to  create  a local 
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modification  of  the  ionosphere  to  create  coupling  into  and  out  of  earth  de- 
tached modes.  These  methods  suffer  from  the  defect  that  the  coupling 
efficiencies  are  small  compared  to  what  one  might  obtain  from  injection 
by  naturally  occurring  ionospheric  tilts.  In  the  first  place,  the  fraction 
of  the  incident  energy  scattered  from  these  ionospheric  irregularities  is 
generally  small;  secondly,  much  of  the  scattered  energy  is  spread  out  into 
undesired  directions. 

7 

High  power,  high  frequency  radio  transmitters  have  been  used  to 
intentionally  modify  the  electron  temperature  and  density  in  the  F region 
of  the  ionosphere.  Under  the  influence  of  a passing  electromagnetic  wave, 
an  oscillatory  motion  is  imparted  to  the  electrons  in  the  plasma  constituting 
the  ionosphere  absorbing  energy  from  the  wave,  and  consequently  heating 
the  electrons  in  a local  region  creating  irregularities  in  the  ionization  and 
electron  temperature.  As  electron  temperature  increases,  the  plasma 
pressure  increases,  and  the  plasma  tends  to  expand.  However,  the  electrons 
are  constrained  to  spiral  about  the  magnetic  field  lines;  thus,  the  plasma 
expands  along  these  lines. 

Consequently,  the  irregularities  in  the  electron  density  and  temperature 
become  aligned  with  the  local  geomagnetic  field,  and  serve  as  scattering 
centers  for  incident  HF  radiation. 

These  scattering  centers  may  be  modeled  by  long  thin  cylinders,  and 
the  incident  radiation  will  be  coherently  scattered  into  a cone  centered  about 
the  local  magnetic  field  line  with  apex  angle  0 equal  to  the  angle  between  the 
directions  of  the  incident  radiation  and  the  local  magnetic  field.  The  lengths 
of  the  scattering  centers  are  presumed  to  be  of  the  order  of  tens  of  kilometers, 
while  the  wavelengths  of  interest  are  from  0.  1 km  (3  MHz)  to  0.  01  km  (30  MHz); 
under  these  conditions,  the  scattered  radiation  is  effectively  confined  to  the 
surface  of  the  scattering  cone,  and  the  far  field  pattern  R of  the  time  averaged 
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radiated  power  is  given  by 


= sin  y 


where  y is  the  angle  between  the  electric  field  vector  of  the  incident 
radiation  and  the  direction  of  the  scattered  radiation.  This  pattern  is 
relative  to  the  forward  scattered  radiation. 


The  geometry  of  the  situation  is  shown  in  Figure  III.  1,  which  displays 
the  geometry  referenced  to  a spherical  polar  coordinate  system  constructed 
with  the  direction  of  the  incident  radiation  as  the  pole;  therefore,  the  2 field 
of  the  incident  radiation  will  lie  in  the  equatorial  plane.  The  "longitude"  is 
measured  from  the  plane  containing  the  incident  radiation  direction  and  the 
local  geomagnetic  field.  The  angle  then,  is  given  by 


cosy  = sin  0 cos  (§  - f 


and  the  far  field  radiation  pattern  is  given  by 


R (®,  $,  §„)  = 1 - sin2 *©  cos2  ($  - $_) 
E E 


Three  special  cases  will  be  considered:  (a)  plane  polarization  with 
the  electric  field  of  the  incident  radiation  in  the  plane  containing  the 
geometric  field  line  and  the  direction  of  incident  radiation  (§  =0); 

(b)  plane  polarizatior  with  the  electric  field  perpendicular  to  this  plane 
($E  = TT/2)’'  (c)  circular  polarization.  The  far  field  radiation  patterns  are 
respectively  given  by 

2 2 2 
R (0,  ?)  = cos  ® + sin  0 sin  $ 
a 

R^  (0,  5)  = 1 - sin4"  © sin2  I 

Rc(©,$)  = 1 - j sin"  0 = | (1  + cos2®) 


It  is  desirable  that  these  results  be  expressed  in  a local  spherical 
coordinate  system  rather  than  the  ray  coordinate  system.  Choosing  the 
magnetic  field  line  as  the  pole  of  the  coordinate  system,  the  incident 
radiation  direction  is  described  by  (0^  = 9,  cp^  = 0)  and  the  locus  of  the 
scattering  cone  is  given  by  (0^  = 9,  cp^  = cp;  cp  = 0,  2tt).  Referring  to  the 
spherical  triangle  KMK1  in  Figure  III.  I,  we  have, 

2 ^ Z 

cos©  = cos  0 + sin  9coscp  = 1 - sin  9 (1  - cosep) 

sin  © sin  $ = sin  9 sincp 


and. 


R (6,  cp) 

a 

Rb(e,  cp) 

R (6,  cp) 

C 


1 z z 

1 -■J  sin  0 cos  0 [3  - 4 cosep  + cos2cpl 

1 z 

1 - — sin  9 [1  - cos  2cp] 

1 sin2  9 [(1  +3  cos2  9)  -4  cos2  9 cosep  - sin2  9 cos2ep] 


3.  9 
3.  10 
3.  11 


These,  then,  are  the  expressions  for  the  far  field  intensity  patterns 

for  the  three  special  polarizations:  Polarization  in  the  plane  defined  by 

the  incident  radiation  and  geomagnetic  field  line  R ; polarization  perpen- 

3. 

dicular  to  this  plane  R ; and  circular  polarization  R . Note  that  the 
elements  on  the  scattering  cone  are  rpecified  by  the  angle  cp  measured 
from  the  direction  of  the  incident  radiation,  and  that  the  radiation 
patterns  are  relative  to  this  direction. 

These  three  intensity  patterns  are  plotted  in  Figures  2 through  7 for 

incident  angles  9 = 15°,  3C P , 45°,  6C f,  75°,  and  90°  respectively.  Pattern 

R reaches  a minimum  at  cp  = 18CP,  where  pattern  R.  reaches  a maximum, 
a D 

The  minimum  of  R,  occurs  at  cp  = 9 CP  ; and  that  of  R varies  from  cp  = 180® 
b c 

when  9 = 45°  to  cp  = 9CP  when  9 = 90° . 
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A 

£ — direction  of  electric  field  vector 

A 

yn  — direction  of  magnetic  field 


Figure  III-l:  Geometry  of  Scattering  Cone 
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FIGURE  HI. 2.  RELATIVE  INTENSITIES  ALONG  SCATTERING  CONE.  THETA 


FIGURE  III. -4.  RELATIVE  INTENSITIES  ALONG  SCATTERING  CONE.  THETA 


FIGURE  UI.S.  RELATIVE  INTENSITIES  ALONG  SCATTERING  CONE.  THETA 


FIGURE  III. 7.  RELATIVE  INTENSITIES  AL0N6  SCATTERING  CONE.  THETA 


The  worst  case  for  R occurs  at  0 = 45°,  and  for  R,  and  R when  0 = 9CP . 

a be 

Although  R^  and  R^  vanish  in  their  respective  worst  case  minima,  R never 

falls  below  0.  25  (6  db  below  the  forward  scattering  direction).  Assuming 

that  the  polarization  is  not  rigorously  controlled,  then  the  pattern  R^  is  most 

relevant,  representing  an  average  state  of  polarization;  if  the  polarization 

may  be  controlled  in  order  to  optimize  the  radiation  intensity  from  any  given 

region  of  the  scattering  cone,  then  one  can  do  better  than  R . We  conclude 

c 

that  the  attenuation  of  the  scattered  signal  due  to  distribution  over  the  scattering 
cone  is  not  a serious  limitation  to  the  feasibility  of  employing  a region  of  field 
aligned  irregularities  to  cause  a change  of  direction  of  incident  radiation.  The 
most  important  consideration  in  the  overall  feasibility  of  this  method  for  in- 
jecting RF  signals  into  long  distance  earth  detached  modes  is  the  fraction 
e^(0)  of  the  incident  radiation  absorbed  and  reradiated  into  the  scattering  cone, 
which  will  in  general  depend  upon  the  angle  of  incidence,  and  may  also  depend 
on  the  state  of  polarization  indicated  by  the  subscript  p.  The  total  scattering 
efficiency  for  injection  into  the  direction  (0,  cp),  then  is  given  by 

T]  (9,cp)  = e (0)*R  (0,cp)/2ttN  (0)  , 3.12 

r r r r 

i r2TT 

where  N (0)  = — 1 R (0 , cp)  dcp  is  a normalization  constant;  for  the  three 
p 2tt  c0  p 

states  of  polarization  considered  above,  the  normalization  constants  are 
given  by, 


3 . 2 

N (0)  = 1- 
a 

- sin  20 

9 

3.  13 

Nb(0)  = 1- 

1 • 2 
2Sm  9 

9 

3.  14 

N (0)  = 1- 
c 

1 . 2 

3 . 2 „ 

3.  15 

- sin  0 - 

— sm  20 

1 b 
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The  derivation  of  the  expressions  and  for  arbitrary  polarization 
is  straightforward,  although  cumbersome.  The  polarization  is  specified  by 
the  angle  $ measured  in  the  plane  perpendicular  to  the  direction  of  the 
incident  radiation  from  the  plane  containing  the  incident  direction  and  the 
focal  magnetic  field  line.  The  expression  for  the  radiation  pattern  becomes: 


R (® > $E) 


1 - sin 


2 


+ 

+ 


cos  $ cos 

E 

2 2 
sin  § sin  $ 

sin  $ cos  $ sin 


3.  16 


A little  trigonometric  manipulation  yields: 


where 


2 2 1 2 2 2 
sin  0 cos  $ = — sin  0 cos  0 f (cp) 


f (p)  = 3-4  cosp  + cos2p 


3.  17 

3.  18 


The  radiation  pattern  then  becomes: 

1 2 2 2 2 
R (0,  p)  = 1 sin  0 [cos  0 f (p).cos  $ 

p L sLj 

+ (1  - cos2p).  sin^  $ + Jz  cos0.  sinp.  f(p).  sin2$  ] 3.19 

Ej  hj 

And  the  normalization  factor  becomes: 

12  ? 2 2 

N (0)  = 1 - — sin  0 [3  cos  0 cos  $ + sin  § ] 3.  20 

p c*  tej  hj 

The  locus  of  the  scattering  cone  may  be  transformed  from  the  local 
geomagnetic  coordinate  system  into  the  local  horizon  coordinate  system. 
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specified  by  elevation  and  azimuth,  in  which  the  criterior  for  earth  detached 
modes  (EDM)  of  propagation  are  defined,  presumably  by  upper  and  lower 
elevation  limits  at  each  local  azimuth.  The  intersections  of  the  scattering 
cone  with  this  band  will  in  general  select  two  windows  into  which  EDM  is 
injected  from  the  incident  radiation.  The  intensity  patterns  may  be  inte- 
grated across  these  windows  to  determine  the  fraction  of  the  incident  radi- 
ation, relative  to  the  forward  scattered  radiation,  injected  into  EDM. 

A problem  not  considered  in  the  present  analysis  is  the  question  of  the 
coupling  efficiency  of  the  field  aligned  irregularities.  This  would  be  ex- 
pressed as  the  fraction  of  the  incident  energy  scattered  into  the  forward 
direction  of  the  scattering  cone. 

The  problem  of  the  criterior  of  injection  and  sustaining  radiation  in 
EDM  propagation  is  considered  in  detail  elsewhere  in  this  report.  This  is 
a difficult,  and  to  some  extent  undefined,  problem  depending  on  the  global 
electron  density  distribution,  the  operating  frequency,  the  (arbitrary) 
minimum  range  qualifying  for  EDM,  the  (arbitrary)  minimum  altitude  of 
ray  perigee  qualifying  for  EDM. 

Clearly,  the  weakest  criterion  is  to  determine  the  elevation  angle  for 
the  straight  line  (high  frequency  limit)  tangent  to  the  earth  from  the  elevated 
source,  and  specify  the  elevation  limits  as  between  plus  and  minus  this  angle. 
These  limits  may  be  further  restricted  by  raising  the  minimum  ray  perigee 
height.  Ultimately,  the  limits  are  determined  by  detailed  ray  tracing  through 
a model  ionosphere.  In  the  study  described  below,  the  criterior  selected  is 
that  radiation  incident  upon  an  elevated  scattering  with  a local  elevation  angle 
up  to  plus  or  minus  5°  is  presumed  to  arrive  via  long  distance  earth  detached 
propagation.  Clearly,  this  presumption  is  subject  to  qualification  based  upon 
analysis  of  the  results  of  the  ray  tracing  program. 
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Given  the  existence  of  a heated  region  of  the  ionosphere,  it  is  of  interest 
to  determine  the  locus  of  transmitter  sites  suitable  to  inject  RF  signals  into 
earth  detached  modes  of  propagation  via  scattering  from  the  field  aligned 
irregularities.  This  problem  may  be  attacked  by  considering  incident  radi- 
ation on  the  heated  regions  from  various  azimuths  at  small  local  elevation 
angles,  selecting  a set  of  scattered  rays  along  the  scattering  cone,  and 
tracing  these  rays  to  the  ground  or  penetration  of  the  ionosphere.  Those 
rays  that  strike  the  ground  trace  out  the  locus  of  suitable  transmitter  sites, 
from  the  principle  of  reciprocity. 


This  ray  tracing  simulation  program  was  carried  out  for  frequencies  of 
6,  12,  18,  24  and  30  MHz  incident  on  field  aligned  irregularities  located  above 
Plattville,  Colorado  at  altitudes  of  150,  200,  and  250  km.  The  radiation  was 
presumed  to  have  come  from  earth  detached  propagation  modes,  incident  from 
magnetic  azimuths  0“ , 45°,  9C f,  135°  andl8Cf,  at  local  elevation  of  Cf  . Ele- 
vations of  ± 5°  were  also  included  at  the  200  km  level. 

The  results  of  this  simulation  are  presented  in  Figures  III.  8 through 
III.  17.  These  charts  trace  the  locus  of  those  rays  on  the  scattering  cone 
which  strike  the  ground  within  the  limits  of  the  charted  region;  sixteen  rays 
were  traced  on  each  scattering  cone.  The  numbers  annotated  to  each  curve 
represent  the  magnetic  azimuth  from  which  the  incident  horizontal  (or  nearly 
horizontal)  ray  comes;  and  therefore,  from  reciprocity,  the  magnetic  azimuth 
into  which  radiation  will  be  horizontally  scattered  from  a transmitter  located 


on  the  charted  curve. 


Figure  III.  8 shows  the  geometry  of  the  scattering  cones,  without  re- 
fraction, in  the  vertical  plane  containing  the  local  geomagnetic  field  line 
appropriate  to  Plattville,  Colorado  where  the  local  geomagnetic  dip  is 
approximately  67° . The  rays  labeled  000M  and  180M  are  incident  horizontally 
from  the  indicated  magnetic  azimuth  on  a scattering  center  at  a height  h above 
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the  surface  of  the  earth;  the  scattering  cones  generated  by  these  rays  are  as 
indicated.  The  ray  from  this  elevated  source,  just  grazing  the  earth  mag- 
netically south  of  Plattville  must  leave  the  source  at  a local  elevation  -8 

Mm 

and  grazes  the  earth  at  a ground  range  R 0 from  Plattville,  where 

o m 

cosg  = R /(R  + h) 

•m  o o 

This  ray  lies  on  the  scattering  cone  generated  from  a horizontally  incident 
ray  from  magnetic  azimuth  z^  given  by 

c°szm  = COS  (A  + 0 )/ cos  A 

For  h«200  km,  0 « 1 4* , and  z «066M. 

m m 

Figures  III.  9»  III.  10,  and  III.  11  show  the  locus  of  the  intersection  on 
the  earth's  surface  of  the  rays  on  the  scattering  cone  after  refraction  within 
the  midnight  ionosphere  in  the  neighborhood  of  Plattville.  In  Figures  III.  9 
and  III.  10,  for  scattering  centers  at  150  km  and  200  km  respectively,  the 
scattering  centers  lie  below  the  bottom  of  the  nighttime  ionosphere  (which, 
of  course,  is  not  technically  possible  for  the  presumed  method  of  creation) 
and  the  downward  scattered  rays  are  unrefracted;  a single  locus  is  generated 
for  all  frequencies.  These  loci  are  the  090M,  135M  and  180M  curves.  At 
6 MHz,  portions  of  the  000M,  045M,  and  090M  scattering  cones  undergo 
refraction  in  the  F layer  and  a subsequent  return  to  the  ground. 

As  the  scattering  center  is  raised  to  250  km  the  effects  of  refraction  in 
the  ionosphere  become  evident  as  shown  in  Figure  HI.  11,  most  clearly  in 
the  090M  traces.  Note  also  that  the  6 MHz  curves  at  azimuths  090M,  135M 
and  1 80M  become  closed,  and  that  the  000M  and  045M  traces  appear.  Noting 
that  the  undeviated  incident  ray  is  a member  of  the  scattering  cone,  a closed 
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curve  clearly  implies  that  this  height  is  a normal  reflection  height;  this 
inference  is  again  reinforced  by  noting  that  the  undeviated  rays  in  the  000M 
and  045M  scattering  cone  also  strike  the  ground  after  refraction  in  the  iono- 
sphere. Thus,  250  km  is  the  normal  reflection  height  for  a 6 MHz  ray 
launched  from  the  ground  returning  to  the  ground  after  reflection  from  the 
ionosphere  at  a range  of  approximately  1420  km  from  the  transmitter. 
Therefore,  we  conclude  that  the  6 MHz  ray  arriving  horizontally  at  the 
250  km  level  is  not  a candidate  for  long  distance  earth  detached  propagation. 
However,  a 6 MHz  signal  so  arriving  may  have  some  energy  scattered  into 
some  appropriate  elevation  angles  for  earth  detached  propagation. 

Figures  III.  12  through  III.  16  show  the  intersections  of  the  scattering 
cones  after  refraction  from  the  Autumn  Noon  ionosphere,  and  Figure  III.  17 
is  the  equivalent  of  Figure  III.  13  for  the  Summer  Noon  ionosphere. 
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